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CURRENT  AND  MAGNETIC  FIELD  DISTRIBUTION 
FOR  AN  INFINITELY  l  ONG  SUPERCONDUCTOR 
OF  RECTANGULAR  CROSS-SECTION 


G.  (Sillier 


I.  INTRODUCTION 

The  nvimerical  results  contained  in  this  report  give  an  essentially 
complete  picture  of  the  quantitative  manner  in  which  current  piles  up 
along  the  corners  of  a  superconductor.  We  have  chosen  various  cross- 
sections  varying  down  to  dimensions  comparable  to  the  skin  depth  as  we 
will  later  wish  to  compare  the  results  obtained  here  (using  the  local 
London  equation)  with  those  to  be  obtained  from  the  non-locad  Pippard 
equation. 
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n.  TWO  REGION  PROBLEM 

We  assume  that  inside  the  superconductor  current  density  and  mag¬ 
netic  field  are  related  by  the  equations 


V  X  (.x'i)  =  H 

Vx  H  =  3 


(London) 

(Maxwell) 


(II.  1) 


It  is  convenient  to  choose  a  rectangular  coordinate  system  with  the  Z  -axis 

■■A 

coinciding  with  the  longitudinal  axis  of  the  superconductor.  Let  A  be  a 

-A 

vector  potential  such  that  H  =V  x  A  with  the  gauge  V  *  A  =  O,  then  the 
only  non-vanishing  component  of  A  is  A^  which  we  will  write  as  the 
scalar  function  ^  (x,  y).  The  current  density  and  the  magnetic  field  are 
given  in  terms  of 


Using  the'  Maxwell  equation  outside  the  superconductor  and  combining 
£qs.  (II.  1),  we  have 


o 

outside 

J-  ^ 

inside, 

(n.3| 


where  the  field  outside  is  still  given  by  (11.2).  At  the  superc  onductor 
boundary  Hx  and  Hy  must  be  continuous;  this  condition  is  fulfilled  with  no 
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loss  of  generality  if  ^  and  its  normal  derivative  are  required  to  be  con¬ 
tinuous  across  the  boundary. 


Since  (U.  3)  is  homogeneous,  the  normalization  of 
at  will.  Now  far  away  from  the  conductor,  IHI  = 
the  total  current  is 


may  be  chosen 
so  from  (II.  2) 


(11.4) 


Our  mathematical  and  computational  problem  is  then  to  solve  the  two- 
regional  problem  (II.  3)  by  a  function  which  is  continuous  and  has  continu¬ 
ous  first  derivatives  everywhere  and  which  is  asymptotically  logarithmic. 
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I 

1 


m.  REGIONAL,  ITERATION 

The  iterative  method  selected  for  solving  (II.  3)  consists  of  alter¬ 
nately  solving  an  exterior  and  Interior  problem.  Let  C  denote  the  boun¬ 
dary  of  the  rectangle  and  suppose  (c-)  is  an  initial  guess  for  the  boundary 
.  Let  ^  ^  and 

terior  problems  respectively. 


values 


ofcp 


denote  solutions  of  the  interior  and  ex- 


The  iteration  is  started  by  solving 

and  continued  by  solving 


(m.i) 


(ni.2) 


Since  we  must  ultimately  match  both  the  values  of  the  functions  and  their 
normal  derivatives,  there  are  many  ways  to  choose  the  boundary  condi¬ 
tions  other  than  those  selected  above.  To  see  that  this  choice  is  not 
arbitrary,  consider  a  one -dimensional  analogue  of  (11.3), 

d  ^ 

=  O  ,  ixl  >  t 
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To  effeet  a  normalization  of  set  ^(KL)  =  1,  K>  l,then 

^CX)  = ,  gpcx)  *  1  4.3(x-ICt>  ^  (m,4, 

According;  to  our  proposed  iterative  method,  we  compute 

Ah.  (V)  =  Ba 

(in.  5) 

I  *  B»,.(l-K)L  =  AftCoJcC-r? 


■D  _  X  c/r^  (-^}p 


and  the  condition  for  convergence  is 

A  QjcrtJi^'hc)^  ✓  1 

L  K-1 

In  our  problem,  we  are  principally  interested  in  large  K  so  this  is  the 
appropriate  order  of  the  iteration.  Reversing  the  above  procedure  leads 
to  the  reciprocal  of  the  above  condition;  this  would  be  the  right  choice 
tpr  K  near  1  . 
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IV.  EXTERIOR  PROBLEM 

Since  ^  is  known  to  be  constant  on  sufficiently  large  circles,  one 
could  conceivably  specify  on  such  a  circle  and  attempt  to  solve  the 

exterior  problem  by  a  relaxation  method  for  some  difference  approxima¬ 
tion  of  Laplaces  equation.  There  are  two  substantial  drawbacks  to  this 
approach;  a)  because  of  the  accuracy  required  near  the  rectangular 
boundary  an  enormous  number  of  grid  points  would  be  needed,  b)  there 
is  no  satisfactory  way  to  determine  how  large  the  bounding  circle  should 
be  (it  turns  out  that  it  would  have  to  be  exceedingly  large  for  long,  thin 
rectangles).  With  the  help  of  the  conformal  mapping  developed  below,  we 
are  able  to  reduce  the  problem  to  series  form  and  thereby  treat  the  whole 
infinite  region. 

Consider  a  quadrant  of  the  exterior  as  shown  below. 


Figure  IV.  1  Quadrant  of  the  Region  Exterior  to  the  Rectangle 
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The  region  A  is  the  conformal  image  of  a  half strip  in  the  w-plane 
under  the  mapping  (c.f .  Appendix  A), 


2  Cw;  -  /V  cxH,  X  -  af 

Where  W  =  NB(o(),  L  =  NA(flt)  and 

Syan,t  -CH^  d'c,  ^ 


On  this  half-strip,  the  solution  of  Laplaces  equation  is  given  by  the  series 


^  (  U  ,  v)  =  ^  ^  ^  ^  c 


-niA 


(IV.  4) 


and  hence  ^  (xiY)  =  i  (a  (x,  y),  V  (x,  y)^  is  obtained  by  inverting  (IV.  2). 

The  coeffieient  of  v  in  (IV. 4)  determines  the  normalization  of  the  asymptotic 
form  of  r  •  On  x  =  o,  we  have 


v-Nfc^-iae^)  *,  ^  j  -  -S: 


(IV.  5) 


and  hence  1  -  4ir  is  the  total  current  for  all  cases  computed. 
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V.  INTERIOR  PROBLEM 

The  solution  for  the  interior  problem  was  approximated  by  solving 
a  difference  equation  using  a  relaxative  method.  A  rectangular  grid  was 
selected  with  a  fine  subdivision  near  the  material  boundary  and  with  a 
comparatively  coarse  subdivision  near  the  middle  of  the  superconductor. 
Error  estimates  were  obtained  for  this  interior  problem  both  from  known 
local  formtilas  and  by  computing  a  test  problem  of  known  solution. 


VI. 


COMPUTATION 
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These  computations  were  performed  on  the  Space  Technology 
Laboratories'  I.  B.  M.  7090  and  required  some  five  to  ten  minute  computa¬ 
tion  time  per  case.  In  practice,  it  was  found  that  only  a  few  regional 
iterations  were  needed  and  that  convergence  was  hastened  by  averaging 
the  coefficients  in  successive  exterior  series  (IV.  4).  For  a  typical  case, 
several  thousand  grid  points  were  used  in  the  interior  region  and  several 
hiuidred  terms  were  used  in  the  exterior  series.  At  this  time  it  appears 
to  the  author  that  a  series  approach  to  the  interior  problem  would  have 
significant  advantage  over  the  relaxation  method  used. 

Results  of  the  computation  are  exhibited  in  the  graphs  of  Appendix  B. 
For  large  regions,  the  graphs  only  show  the  character  of  ^  near  a  corner. 
The  labeling  fc-  fn-  refer  to  the  values  ^  (L+  iW), 

^  (L)  as  shown  in  Figure  (IV.  1). 
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APPENDIX  A 

The  mapping  used  in  the  exterior  problem . 

Let  Z(W)  be  the  mapping  function  for  the  regions  shown  below. 


W  -  PLANE 


Z  -  PLANE 


Figure  A.  1  Regions  Defining  Z  (W) 


We  will  show  that 


Z  (w)  =  i  B  +  5'v^T' 

O' 

A  =  <A-2> 

The  fundamental  net  for  the  region  exterior  to  a  rectangle  consists  of  curves 
y  (x,  y)  =  ctmstantyv  (x,  y)  =  constant,  where  u  and  v  are  coi^ugate 
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harmonic  functions  on  this  region.  Cutting  this  region  along  the  lines  of 
symmetry  results  in  the  region  Jd!  whose  boundary  ie  entirely  made  up 
by  lines  of  the  fundamental  net*  thus  the  mapping  W  (Z)  =  u(x,  y)  -f  iv(x,  y) 
will  map  JCV'  on  to  a  region  A  with  a  rectangular  fundamental  net  and  we 

can  determine  the  image  SL  o£  SL*  under  W  (Z)  as  well  as  JTl"  under 

dW  • 

^  ^^x  "  ^y  “  ^  **i*PP^ 

F(W),  we  will  obtain  a  differential  equation 

i?  ■  ■  H-  =  ifci 

which  can  be  integrated  for  the  desired  mapping  2  (W). 

The  region  Jfl!*  is  constructed  by  considering  the  values  of  u^^t 
Uy,  Vy  on  the  various  pieces  of  the  boundary  of  JTI  . 


Figure  A>  4  The  Image  of  flf  under 
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The  image  SI 


Ml 


M  2 

of  JTL  under  the  mapping  l/L  is.: 


M"  (A)”  (A+iB)"  (iB)"  (00)" 


Figure  A.  5  An  Intermediate  Image 

Now  choose  Gj  =  it,  ~  O,  =  O,  and  map  by  cob^w)  on  to  St  , 
this  gives 


-I  I 


Figure  A,  6  An  Intermediate  Image 

which  except  for  translation  and  normalization  agrees  with  -O-  *  The 

normalization  is  fixed  by  (iB)*'  -  (A)*'  =  2  and  since  (A  +  iB)"  lies  between 
(A)"  and  (iB)"  we  can  choose  the  translation  to  be  cos  for  0  <  o(  <  ir, 
then 


1- 

L"* 


Gm,  W  -  Co^  oc 


(A.7) 
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But  1  /  L  = 


so  Z(w)  is  determined  by 


diE 

dw 


(A.  8) 


This  can  be  developed  in  an  exponential  series  by  using  the  generating  ftinc- 
tion  for  Liegendre  polynomials 


(A.  9) 


To  do  this,  write 

di  = 

dw 


tW 

JL  A  ^ 

VjT  ^ 


icy 

X  S  ^ 


f- 


(A.  10) 


then 


ziW 

- 


1  (ne'*' to)  =  ^  -$ - 

dv»/ \  6^wl  yx- zc^oCe^-i- e 

_  .  (  ziw  iU/N  ^  c'ytiy^  ,  V 


(A,  11) 


Iff* 


=i  2  (zi^>  -  \ 

i  #•* 


This  gives 
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aw  vr^  /r  ®  (A. 


12) 


4  «  -  _ 

but  f  or  e  ^  near  aero,  ^  C  =  1, 
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APPENDIX  B 

The  following  graphs  exhibit  the  lines  of  constant  current  density 
in  the  superconductor  and  lines  of  constant  magnetic  field  outside  the 
superconductor. 
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